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Various properties of mesoscopic two-dimensional Josephson junction arrays are reviewed. Particular attention 
is paid to structure of the topological excitations, charges and vortices, which are shown to be dual to each other. 
This duality persists in the presence of external magnetic fields and offset charges, which influence vortices and 
charges in an equivalent way. A double-layer junction array is also considered, where an even further reaching 
duality is discovered. 



1. Introduction 

Since the original work by Kramers and Wan- 
nier Q on the two-dimensional (2D) Ising model, 
duality has been proven to be a powerful tool in 
field theory and statistical mechanics 0]. The 
idea behind this transformation is the mapping 
of the weak coupling region of the system under 
consideration onto the strong coupling range (and 
vice versa). The symmetries of the system un- 
der this transformation lead to important insight 
into the structure of the model, especially in the 
intermediate region of coupling constants which 
is usually elusive to standard treatments. Dual 
transformations applied to the topological excita- 
tions of the system are particularly advantageous 
since it is possible to recast the partition func- 
tion solely in terms of these degrees of freedom. 
Jose et al. M applied these techniques to show 
that the vortices in 2D XY-model can be mapped 
onto the charges of a two-dimensional Coulomb 
gas. Kadanoff m showed that this mapping is not 
restricted to the XY-model, but various systems 
can be mapped onto coupled Coulomb gases. A 
review of these techniques applied to a number of 
systems can be found in Ref. ||. 

In the past years two-dimensional Josephson 
junction arrays (JJA) have proven to be an ex- 
cellent arena for the study of a variety of phase 
transitions J^] . They are fabricated from an array 
of superconducting islands connected by Joseph- 



son links. Each island is characterized by the 
modulus and the phase of the order parameter, 
Aexp(i0). Upon lowering the temperature, each 
island of the array goes superconducting at the 
BCS critical temperature T c0 . But, in spite of 
the fact that each island is superconducting, the 
whole array remains in a resistive state as long as 
the phases, driven both by thermal and quantum 
fluctuations, have not acquired (sufficient) long 
range order. This global phase coherence sets in 
at a lower temperature. It is reasonable to assume 
that at these temperatures the magnitude of the 
order parameter is fixed to its equilibrium value, 
and the array can be described only in terms of 
the phases <pi. The resulting Hamiltonian is 

H=-E J J2 cos ( ( l>i-<f > j) . (!) 

m 

where Ej > is the Josephson coupling energy, 
and the summation is over nearest neighbors. 
The lattice model defined in Eq.(Q) is equivalent 
to the 2D XY-model - a planar lattice of localized 
planar magnetic momenta Si of unit length. It de- 
scribes other physical systems as well, including 
planar ferromagnets, two-dimensional crystals |t]] 
or two-dimensional superconducting films B. 

The XY-model undergoes a Berezinskii- 
Kosterlitz-Thouless (BKT) transition 
driven by fluctuations in the vortex density. For 
the JJA this implies that the array is phase coher- 
ent (superconducting) below the temperature Tj, 
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which is of the order of Ej (we use units where 
h = fcg = 1). This is demonstrated by a sim- 
ple argument due to Kosterlitz and Thouless. At 
zero temperature the lowest energy state of the 
system has the phases 4>i aligned, however, at fi- 
nite temperatures the free energy F = E — TS 
should be minimized. Consider an isolated vor- 
tex configuration. Its energy is readily calculated 
from Eq. (Q). Logarithmic divergences are cut 
off at short distances by the lattice spacing a and 
at long distances by the system size R, hence 
E = ttEj \n(R/a). The entropy of the vortex con- 
figuration is given by the logarithm of the number 
of possible positions of the vortex, S = 2 \n(R/a). 
Thus, the free energy of the vortex becomes neg- 
ative at the temperature Tj = irEj/2. Below this 
temperature the phases are ordered, and vortices 
may appear only in bound pairs. I.e. the XY- 
model is in a ferromagnetic state, while the JJA 
is superconducting. Above the transition temper- 
ature Tj free vortices may form. Their motion 
causes dissipation, and hence the state is resis- 
tive. This analysis is confirmed by the renormal- 
ization group calculations fll| . A thorough re- 
view concerning the statics and dynamics of the 
BKT transition can be found in Ref. ]T^ ]. 

2. Charge- vortex duality in 2D Josephson 
junction arrays 

If the superconducting islands are of sub- 
micron size, as can be realized with modern nano- 
lithography, the electrostatic charging energy as- 
sociated with non-neutral configuration of the is- 
lands cannot be disregarded any longer. It can 
be expressed as 

ij 

Here the summation is over all islands i with lo- 
cal voltage Vi = dt(4>i/2e). The capacitance ma- 
trix is determined in a good approximation by 
the on-site and the nearest neighbors elements, 
Cu — Cq + AC, Cij — —C if i and j are nearest- 
neighbors, and C^ = otherwise. Here Co is the 
capacitance of the grains to the ground while C 
is the junction capacitance. 



In a quantum mechanical treatment, the charge 
Qi on the island i has to be viewed as conjugated 
to the phase 4>i ( see e -g- Ref. i- e - 

1 d 

Qi = ld(^/2eY (3) 
Then the Hamiltonian of the array is |l4| 

»3 (ij) 

The interaction between the charges in the array 
is determined by the inverse capacitance matrix 
C^j . Its range depends on the ratio C/Co; if the 
capacitance to ground is much larger than the 
junction capacitance the charges feel an on-site 
interaction. In the opposite limit the Coulomb 
interaction depends logarithmically on the dis- 
tance. The capacitance to the ground Cq leads 
to a screening on the length scale ~ a^/C/Co- 

The two characteristic energy scales in the sys- 
tem are the Josephson coupling energy Ej , which 
is associated to the tunneling of Cooper pairs be- 
tween neighboring islands, and the charging en- 
ergy Eq — e 2 /(2max{Co,C}), which is the en- 
ergy cost to add an extra Cooper pair on a neutral 
island. The charging energy Eq tends to inhibit 
Josephson tunneling and is responsible for quan- 
tum fluctuations of the phases 4>i on each island. 
If Ej ^ Ec the fluctuations of the phases are 
weak, and the system acquires global supercon- 
ducting coherence. 

In this section we derive some properties of 
quantum JJA from duality arguments. Our start- 
ing point is the partition function, 

z= nr^ o) e 

i J " {m;=0 : ±l,...} 
rc/> i0) +2irmi 

/ 1 £>&(r)exp(-S{0}). (5) 

Here the Euclidean path integration is carried out 
with the boundary conditions 

4n(o) = 4 0) ; MP) = 4° ] + ^m t , (6) 

with [3 being the inverse temperature. These non- 
trivial boundary conditions express the fact that 
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the charges of the grains are integer multiples of 
2e The Euclidean effective action S{(j)} 

has the form 



dr 



C 



-tfj^cosfa-fc) L (7) 

«> j 

The model has been studied extensively to 
derive the phase diagram by means of mean 
field H-Jol, Monte Carlo |j, variational || 
and RG [23] calculations. Non-perturbative fea- 
tures of this model and its connections to planar 
Chern-Simons gauge theories were discussed in 
Ref. HU]. Quantum fluctuations lower the tem- 
perature of the BKT phase transition separating 
resistive and superconductive phases ||| below 
the transition temperature Tjo of the classical ar- 
ray (with Eq =0). In the case Co = the shift 
of the transition temperature is 
„2 



T '' = Tj ° " 48^ (8) 
Beyond a critical value of the charging energy Eq 
the transition temperature vanishes, and the ar- 
ray remains insulating even at zero temperature. 

The purpose of this Chapter is to show that 
there is a dual transformation relating the clas- 
sical limit, Ej 3> Eq, to the opposite quantum 
limit, Ej < E G In the latter case 

the quantum fluctuations of the phases are strong, 
and vortices are ill-defined objects. However, in 
this regime the charges on the islands are well- 
defined variables. In the extreme limit Co <C C, 
the interaction between the charges is logarith- 
mic, in the same way as that of the vortices in 
the classical array. In this case the charges form 
a 2D Coulomb gas and undergo a BKT transition 
at temperature Tq ~ Eq, where the state below 
the transition is insulating. If the capacitance to 
the ground is larger, no finite-temperature phase 
transition arises. However, we can still expect a 
crossover from a low-temperature phase with ex- 
ponentially low conductance to a resistive phase 
at finite temperature. 

In order to describe an intermediate situation 
we need a formulation in terms of both charges 



and vortices. Following Ref. Qj27| we first in- 
troduce the island charges in a path integral rep- 
resentation. In terms of phase trajectories 4>i(t) 
and charges qi(r) = Qi(r)/2e the partition func- 
tion takes the form 

z = Uf ^ I D qj (r)l[ r d<j>f J2 

D^(r)exp(-5{^,g}), (9) 

where the phases obey the boundary conditions 
(^), while the charge paths are periodic, qj(0) = 
qj(/3) — qjo. The effective action in terms of 
phases 4>i( T ) an d charges qi(r) is 

S{hq} = J 'dr j 2e 2 5>C£ 1 ?J - 

+ i ^2 <n4>i - e j E cos (^ ~ fa) \ ■ ( 10 ) 

(ij) j 

The summation over winding numbers {m,} fixes 
the charges to be integer- valued flEf . 

Following the steps discussed in the Appendix, 
the partition function can be expressed as a sum 
over charge and vortex configurations 



Z = E E cx P (-S{g,«}) . 

{q*r}{v,-r} 



(11) 



The effective action for the integer charges qi(r) 
and vorticities Vi(r) is 

r0 , 

S{q,v} = drV 2e 2 q i C~ j 1 q ] + wEjViGijVj 

+ iqSijVj + 4^Bj? iG! y*}- ( 12 ) 

It describes two coupled Coulomb gases. The 
charges interact via the inverse capacitance ma- 
trix. The interaction among the vortices is de- 
scribed by the kernel Gjj, which is obtained as 
the Fourier transform of k~ 2 . At large distances 
i\j 3> a between the sites i and j it depends log- 
arithmically on the distance 



(13) 
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The charges and vortices are coupled in the dy- 
namical theory by the third term. Here 



Qij = arctan 



- Uo 



(14) 



describes the phase configuration at the site i if 
a vortex is placed at the site j. The coupling has 
a simple physical interpretation: a change of vor- 
ticity at site j produces a voltage at site i which 
is felt by the charge at this location. The last 
term qGq represents a spin-wave contribution to 
the charge correlation function. 

The effective action ( |l2| ) shows a high degree of 
symmetry between vortex and charge degrees of 
freedom. In particular, in the limit Cq <C C the 
inverse capacitance matrix depends on distance 
in the same way as the vortex interaction, 



,2 r -l _ E C r 



(15) 



and charges and vortices are (nearly) dual. The 
duality is broken by the last term qGq. This term 
is "irrelevant" for the phase transitions, i.e. it 
merely shifts the transition point. But it has im- 
portant consequences for the dynamical behavior. 

The action ( |l2] ) is written in terms of integer 
charges and vorticities. However, depending on 
the coupling constants, only one type of excita- 
tions may be well defined. In the quasiclassical 
limit, Ej Eq, vortices are the relevant excita- 
tions, while charges are strongly fluctuating and 
can be treated as continuous variables. By inte- 
grating out the charges, one obtains the effective 
action for vortices Vi(r) 

( (i ( 1 

ijkl 

■ ;:/•;., ^r//,,,-,}. (16) 

ij 

In the limit of low self-capacitances, Co <C C, the 
kernel in the first term becomes QjkG^Qij, and 
the effective action for vortices reduces to |2g] 



nEj^ViGijVj} . (17) 



The summation in the partition function is con- 
strained by the neutrality condition ^\ Vi — 0. 

The effective action ( [1^ ) describes a quantum 
Coulomb gas of vortices. This becomes clearer if 
we change from a description in terms of the vor- 
ticity at site i to a continuous description where 
we label the vortices by their center coordinate 
r(r) and sign of vorticity v n = ±1 Both are 
related by 

n 

In these variable the partition function becomes, 

oo . 

Z = E r ^(r) exp(-S{r}) . (19) 

(We have explicitly used the charge neutrality 
conditions: the integration is carried over N vor- 
tices with v = 1 and N anti-vortices with v = 
— 1). The effective action expressed in terms of 
the vortex coordinates is 



S{r} 



dr 



2N 

E 

m,n— 1 



+ ttEjG( 

I'm *V) 



(20) 



The second term in ( pOf ) is the interaction energy 
of the vortices. In the first term we introduced 
the vortex mass tensor |3(J] , 



M a p{r) 



4£ c 



V a V^G(r). 



(21) 



For r 3> a it decreases as r 2 , and consequently 
may be approximated by a local function 



M Q/3 (r) = M v 5 a0 5(r), M v = 



Aa 2 E, 



c 



(22) 



It defines the vortex mass M v [|3JJ,|3l|. In this 
case the first term represents a kinetic energy of 
vortices, which can be rewritten as 



^fi = / d - E ^(r) 

Jo m=l 



(23) 
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A similar consideration for Eq 3> Ej lead us 
to the action of a 2D Coulomb gas of charges 
with charge mass M q = (a 2 E])~ 1 (see also Ref. 
p2[). Thus, in two dimensional arrays of Joseph- 
son junctions a charge-vortex duality exists. In 
the limit Ej 3> Eq the vortices are well-defined. 
They form a Coulomb gas, and can be consid- 
ered as particles with masses. In the opposite 
limit Eq <C Ej the charges are the relevant exci- 
tations. The charges have the same properties as 
the vortices in the corresponding limiting cases. 

3. Double-layered arrays 

Another system of increasing experimental in- 
terest is composed of two parallel 2D Josephson 
junction arrays with purely capacitive coupling 
between them (no Josephson coupling (3j|). We 
will show that in this case an even further reach- 
ing duality between charges and vortices arises. 
We restrict ourselves to the most interesting sit- 
uation, when one array is in the quasi-classical 
(vortex) regime while the other is in the quantum 
(charge) regime. Then the vortices in one layer 
and the charges in the other are simultaneously 
well-defined dynamical variables. (In contrast in 
a single array vortices or charges acquire a kinetic 
energy only after the other variable is integrated 
out.) Another important feature of the present 
system is that the strength of interaction between 
charges and vortices is controlled by the interlayer 
coupling C x and consequently may be tuned inde- 
pendently. We show that the physical realization 
of this interaction is rather different from that in 
a single array. The theoretical description of this 
system has been developed in Refs. 
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The partition function of the system in terms 
of the phases 0^ (the indices i label the islands 
in each array and ^ = 1,2 refers to the layer) is 



2ji 



{m il },{mi2} 



J D^ 1 (r) J D</» l2 (r)exp(-5{0 1 ,0 2 }) , (24) 



with boundary conditions 0^(0) = (pff} and 



(ft) 



%^1 

bf^ L + limiip. The Euclidean effective 



action S{4>i, $2} has the form 



H=l,2 



2 - E Jfl ^2 C0S (<?V - fan) 



8e 2 



£( 



;) 2 



(25) 



Here Co^ are the capacitances of the islands in the 
array fi relative to the ground, are the capaci- 
tances of the junctions in the array /i, and C x are 
the interlayer capacitances between adjacent is- 
lands, while Ejn are the Josephson coupling con- 
stants in the layers. 

We concentrate on the situation in which the 
array 1 is in the quantum (charge) regime while 
the array 2 is in the quasi-classical (vortex) 
regime, i.e. 

E n < e 2 /2C u E i2 » e 2 /2C 2 , 

with C M = maxjCop, C M , C x }. Various types of 
problems van be studied in these systems. If we 
are interested only in behavior of the vortex ar- 
ray (e.g. in an experimental realization where one 
measures the quantities in the vortex array only) 
we can integrate out all degrees of freedom re- 
lated to the charge array and study an effective 
action for the vortex array. An analysis analo- 
gous to that presented for a single array shows 
that this action is essentially that of the Coulomb 
gas [{34| . Consequently the vortex array under- 
goes a BKT transition, and in the quasiclassical 
regime its temperature is lowered due to electro- 
static coupling to the charge array. 

A similar problem can be solved for the charge 
array. Then the charge-BKT transition changes 
into a crossover, since the capacitive coupling to 
the vortex array leads to screening similar as ca- 
pacitances to the ground. 

Our purpose here is to consider both arrays si- 
multaneously and to investigate the charge- vortex 
duality in this system. Similar as in a single array, 
we move from a description in terms of phases 
to one in terms of charges and vortices, and 



6 



use the duality of the resulting action to inves- 
tigate the transition. We will show that charges 
and vortices in this system can be considered as 
two-dimensional dynamical particles with masses. 
The charge-charge and vortex-vortex interaction 
are essentially those of 2D Coulomb particles, 
while the charge-vortex interaction is highly in- 
teresting. 

Before we proceed with the calculation it is nec- 
essary to stress the following. In the regime of 
interest the intcrlayer capacitances C x not only 
couple the layers, but also renormalize the ca- 
pacitances Cbi and C02 of the islands to the 
ground. The physical reason for this is that due 
to the strong fluctuations of charges in layer 2 
and vorticities in layer 1 these variables are ef- 
fectively continuous, and hence a coupling to the 
other array plays the same role as a coupling to 
the ground. Hence the interaction between the 
charges in each layer has a finite range for any 
non-zero C x due to the screening, and the BKT 
transition is replaced by a crossover. On the other 
hand, in the limit C01 C C x C Ci the screening 
length £1 ~ a(Ci/C x ) 1 / 2 can be very large. Be- 
low we assume that these inequalities are satisfied 
and the range of interaction £1 is large enough to 
make it meaningful to speak about the charge- 
unbinding transition. This just means that the 
crossover is (exponentially) sharp. For not so 
weak coupling C x this description becomes mean- 
ingless, since the crossover is strongly smeared, 
and the insulating phase is absent. 

To proceed it is convenient to introduce the 
large capacitance matrix 



C = 



C\ — C x 
— C x C2 



(26) 



Here C M is the capacitance matrix in the array 
fi while (Cx) y = C^Sij. The inverse matrix, de- 
scribing the interaction of charges, in the Fourier 
representation reads as 

c-^k) = j n ( c °» + c ^ 2 + a ) - c *) 

l»=l,2 J 



1 C02 + C*2k + C x C x , 



CL 



The effective action (f25|) can be rewritten in terms 
of integer charges q^ (r) and phases (fr^ (r) of each 
island 



S{q,cj ) } = jf dr^2e 2 J^ Qi 



1 fi 



jm 2^ cos( - 



) -(28) 



It is possible to introduce vortex degrees of free- 
dom in the same way as for one array. Since this 
procedure deals only with the phase variables and 
does not affect the charge interaction (the first 
term in Eq. (ps|)), the generalization for double- 
layered arrays is trivial. We obtain 



Z 



£ £ exp(-S{q,v}),(29) 

{9il(T")>9i2(r)} {Oil(T),»fl(T)} 



where the effective action for integer charges q. 



1 it 



and vorticities is 



S{q,v} = ^ dr|2e 2 £ ^(C^ 1 )^ 



■E 



+ TTEj fl F(e IJ ,E JfJ ,) ^ v^G^v^ 



i qifj.QijVjfj, 



(30) 



Here we wrote terms which arise from the dis- 
cretization of the time with lattice spacing (in ar- 
ray fi) of order e M ~ (SE^Ecp,)" 1 ^ and E Cl i = 
e 2 /2C fJl . (See the Appendix for the function F 
renormalizing the Josephson coupling). In the 
limit considered only the Josephson coupling of 
array 1 is noticeably renormalizcd, 



E n ~ {8E. n E cl ) 1/2 [m(£ C i/£ji)] 



(31) 



We suppress the tilde from now on. Later we 
will assume that the linear size of each array 
is much less than the range of interaction £ M = 
a(C / _ I /C x ) 1 / 2 . This means, in particular, that we 
assume C x <C C%. 
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The action ( p0| ) depends on the charges and 
vorticities in both layers. However, in our situa- 
tion, when the layers 1 and 2 are in the charge and 
vortex regimes, respectively, the vortices in the 
layer 1 and the charges in the layer 2 are strongly 
fluctuating degrees of freedom and may be inte- 
grated out. To do this we suppose the latter vari- 
ables to be continuous and neglect the spin-wave 
charge coupling term, q 2 Gq2, in layer 2. Then 
after performing the Gaussian integration we ob- 
tain the effective action for charges qn{r) in layer 
1 and vorticities ^(t) in layer 2 (to be referred 
below as qi and v{) 



S{q,v} 



dr 



„2 \ " 



8Ee- 



ijkl 



(c- 1 



22 



A / 



•<E*< e * (C- 1 ) u (C- 1 ) ,.1S -(32) 
rrf, L iki i J h 1 

IJKL 



In the limit where the charge interaction is long- 
range we can use the appropriate limits of the 
large capacitance matrix (27). In this case the 
action reduces to 



S{q,v} 



f 2E a ^ 
dr < } qiGijq-j 

o L ^ 



8Eq2 



Gi 



a 



47r 2 CiC 2 



T~7T E QikGklQl 



fei 



+ 



2ttCi 



E ViQikGkjqj 



ijk 



(33) 



This form displays the duality between charges 
and vortices in the appropriate limit. 

The action ( |33| ) is the central result of this 
section. It looks rather similar to the effective 



charge-vortex action in one Josephson junction 
array. The most important difference is that 
in one layer either charges or vortices are well- 
defined degrees of freedom, while the action ( |3^ ) 
of the double-layer array describes the system of 
two well-defined dynamic variables on each site - 
the charges in layer 1 and the vortices in layer 2. 
The action shows a duality between charges and 
vortices (the second term in the square brack- 
ets is small for C x <C 61,62) . Both kinetic 
terms for charges and vortices violate the dual- 
ity due to the numerical coefficients. However, 
close enough to the transitions these terms pro- 
duce only small renormalization of the transition 
temperature, and are irrelevant. Another inter- 
esting feature of this action is that the last term, 
describing the interaction between charges and 
vortices, is also small, while in a single-layer ar- 
ray the interaction is always of the same order of 
magnitude as the other terms. 

It is obvious that for long-range interaction of 
the charges in layer 1 these also exhibit a BKT 
transition, and under the conditions where the ac- 
tion (|33] ) was obtained the transition temperature 
does not feel the presence of layer 2. Hence 



T r =T ( 



CO 



En 
24tt 



To understand the physics described by the ac- 
tion ([33j) it is instructive to map this model onto 
the 2D Coulomb gas. For this purpose we move 
again from the space-time lattice to the continu- 
ous medium and introduce the coordinates of the 
vortex centers and charges 



h(t) -> E 9m ^ r - r m (r)) 

m 

v i{ T ) E^™ S ( T ~ R "( T )) 



(34) 



Here q m = ±1 and v n = ±1 represent charges 
and vortices, respectively, and r m (r) and R„(t) 
the coordinates of their centers. In terms of these 
variables the partition function reads 

oc 00 „ 

Z= E E / Dv 1 {t)...Dv 2M {t) 

M=0 JV=0 

DRi(r) . . . DR 2N (t) exp(-S*{r, R}), (35) 
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with an effective action S{r, R} describing a neu- 
tral system of 2M positive and negative charges 
(q = ±f ) and of 2N positive and negative vortices 
(v = ±l). 

The first and third terms of the action (|3^) can 
be easily transformed by means of decomposition 
(|34|). The first one produces the potential energy 
of charge interaction, 



o(q) _ 

°int ~~ 



2E, 



CI 



dr 



2M 

E 

m,n—l 



q m q n G(r m (T) - r„(r)). (36) 



In principle, the summation includes the terms 
with m = n; these, however, give rise only to the 
chemical potential for charges. The third term in 
Eq. (33) yields the interaction of vortices, 



5; 



(v) 



2 A 



nE J2 / dr 22 v m v n G(R m ( T ) -R„(r)). (37) 

•*° m,n=l 

Here again the term with m — n gives rise to the 
chemical potential for vortices. The terms ( |36] ) 
and ( |37| ) are essentially the action for (classical) 
Coulomb gases of charges and vortices, respec- 
tively @. 

If we neglect the small correction proportional 
to the Cl/CiC 2 in the fourth term in Eq.@ 
then the second and fourth terms can be trans- 
formed to the kinetic energy of charges and vor- 
tices respectively |l4|]. The second term gives 

9mQn^M 7<5 (r m -r n )f*,(38) 



dr 



2M 

E 

m,n— 1 



where the charge mass tensor can be approxi- 
mated as 

M lS (r) = _-jL-V 7 V«G(r) « M q 5 lS S(r) 



1 



, r (39) 

Then the kinetic term for charges takes a simple 
form 



o(q) _ 

D kin — 
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2M 
m— 1 



(40) 



Similarly, the fourth term in Eq. ( p3| ) produces the 
kinetic term for vortices 



o(v) 
D kin 



with 



f P 2N 
•'° m=l 



(r) 



M v = 



Aa 2 E, 



C2 



(41) 



Finally, the last term in Eq.(p3[)) is responsible 
for the interaction between charges and vortices. 
The corresponding term in the action is 



dr 



I 



dr'V R 6(R m -r')G(r'-r„)R m (r) . (42) 



After the integration over r' it reduces to 



cLt Y iv m R m (T)A(R m ) 



with 

A(R m ) = ^g„a(R m (r) - r„(r)) , 



and 



i a 

"Ci 



i 



2 In ('2 



z x r 



(43) 



(44) 



(45) 



The charges are the sources of a vector potential 
felt by the moving vortices. It depends of the 
signs of the corresponding vortices and charges. 
We can rewrite this term, after a partial integra- 
tion in Eq. (|33|) , in such a form that the vortices 
create a gauge potential for the charges. Hence 
this charge- vortex interaction term fully preserves 
the duality between charges and vortices. 
In summary the resulting action 



5{r, R} = S<$ + Sg> + S™ + + S qv (46) 



j(v) 



is essentially that of two coupled 2D Coulomb 
gases. It is symmetric with respect to charges 
and vortices. Both can be considered as particles 
with masses 



1 



M q = 

a 2 E n 



and M v = 



4a 2 £ C2 



(47) 
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respectively. Charges interact via the effective ca- 
pacitance, vortices via the usual logarithmic in- 
teraction with strength Ej 2 - In addition, there is 
a specific feature of two-layered system: the vor- 
tices produce a vector potential a for the charges, 
whereas charges create a gauge potential for vor- 
tices. The magnetic field, associated with this 
vector potential is 

£ = ±-i^ln-, a<r<&. (48) 
2ea z Ci r 

The interaction described by this vector/gauge 
potential is always controlled by the small pa- 
rameter Cyc/C\. In the regime considered it is 
weak. 

4. Magnetic field offset charge duality 

In this last section we return to the case of a 
single Josephson junction array. We assume that 
the array is in a perpendicular position-dependent 
magnetic field Bj. Furthermore, we allow for "off- 
set" charges on the grains Q X i, which in princi- 
ple can be controlled by applied gate voltages be- 
tween the ground and the islands. (In practice 
they are caused e.g. by random impurity charges 
in the substrate.) Both these factors can be ac- 
counted in the formalism described above. The 
action describing this system is (see e.g. |3q| ) 

S{<j>,q} = 



J dr| 2e 2 5>(T) + fc ]Cr 1 feM 



+ i J2 h{t-)Mt) - E 3 cos [^ - fa ~ A v\ \ ■ ( 49 ) 



Here q X i — Q X i/2e, and the magnetic field is 
taken into account by the vector-potential Ay- = 
(2e/c) J? A • dl. After the same transformations 
as those performed in Sec. 2, we arrive to the 
coupled-Coulomb gas action for charges and vor- 
ticities, similar to Eq. ( p"2| ) 

S {q,v} = 

J drj^l 2e ' [* ( r ) + to] C T/kj W + Qxj] ■ 



+nEj[ Vi (T) + [^(t) + /,•] 
+i[q l (r) + q xi \@ijVj(T) + ^^ii( T ) G ij j • ( 50 ) 



The "frustration" 

fi = (2tt) 



describes the magnetic flux through the plaque- 
tte i, measured in units of the flux quantum 
$o = 7rc/2e. Note that due to the fact that the 
charges and vorticities are integer- valued, only 
the fractional parts of both offset charges g X i and 
frustrations fi matters. 

The presence of offset charges Q X i = 2eq X i, 
breaking the particle-hole symmetry, has pro- 
found consequence on the vortex motion. Offset 
charges at positions fi are responsible for a gauge 
potential 



A(r) 



E 



(51) 



which acts on the vortices in the same way as an 
ordinary vector potential acts on charges. Asso- 
ciated with this 'vector potential' A(r) is a 'mag- 
netic field' and a 'Lorentz force' acting on moving 
vortices. We call this force a Magnus force. A 

vortex with vorticity v n and velocity R feels the 
Magnus force fjfl 



F M = v n q x z x R 



(52) 



Here we assumed for simplicity a homogeneous 
gate charge. As a result of the combined effect 
of the Magnus force ( |52|) and the Lorentz force, 
which is the force on the vortex due to an exter- 
nal feeding current, the vortices will move at a 
certain angle with respect to the current. This 
angle is called Hall angle. Its measurement yields 
information on the different forces in the system. 

Real samples are usually characterized by ran- 
dom offset charges. As a result the Magnus force 
averages to approximately zero. This effect is 
probably responsible for the small size of the ob- 
served Hall angle in Josephson junction arrays. 
The forces on vortices in Josephson junction ar- 
rays have been discussed recently in the litera- 
ture, in part in a controversial way (see Ref. p^ - 
fnj). Here we want to stress that only the offset 
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charges, which are responsible for a local devia- 
tions from charge neutrality in the array, lead to 
the Magnus force. 



5. Conclusions 

We have considered two-dimensional Josephson 
junction arrays. In the classical limit, Ej 3> Eq, 
the vortices are the relevant excitations. In the 
opposite quantum limit Eq 3> Ej the charges 
of the islands are well-defined degrees of free- 
dom. Both can be described simultaneously by a 
coupled-Coulomb-gas action, which is dual under 
the interchange of charges and vortices. In each 
of the two limiting cases (Ej <C / S> Eq) the sys- 
tem reduces to a one-component Coulomb gas, 
where either the charges or the vortices can be 
considered as (logarithmically) interacting mas- 
sive quantum particles. Furthermore, the exter- 
nal magnetic field plays the same role for charges, 
as the external ("offset") charges for vortices. 

The situation is even more interesting in a sys- 
tem of two parallel, capacitively coupled Joseph- 
son junction arrays. It is described by a coupled- 
Coulomb-gas action for charges and vortices in 
both arrays. If one array is in the semi-classical 
(vortex) regime, while the other is in the quantum 
(charge) regime we can integrate out the strongly 
fluctuating variables. In this case we arrive at 
an effective action, describing dual charges in one 
array and vortices in the other, which both are 
now dynamic degrees of freedom, in contrast to 
the one-layer problem. Furthermore, vortices and 
charges interact via gauge field of strength con- 
trolled by the interlayer capacitance. This pe- 
culiar interaction between charges and vortices 
resembles the composite fermion scenario of the 
fractional quantum Hall effect, which may be- 
come a subject of future research. 
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Appendix 

In this Appendix we provide the steps leading 
from Eqs. ([§, @ to Eqs. ([□]), (|l|). Starting 
from the partition function (^), we first introduce 
the vortex degrees of freedom. This can be done 
by means of the Villain transformation jl2| (see 
also H); the time-dependent quantum problem 
requires some additional steps [lj,|2^]. We intro- 
duce the lattice with spacing e in time direction; 
this spacing is of order of inverse Josephson fre- 
quency: e ~ {%E 3 E C ) 1/2 - In the Villain approxi- 
mation one replaces 



exp \ - eE, 



[1 - COS(<^ >r - </>,>)]} 



cx p{- 

{m 1T } 



eEjF(eEj) 



J2 |V^ T -27rm IT | 2 }. 



(Al) 

Here we have introduced a two-dimensional vec- 
tor field nii r , defined on dual lattice (alterna- 
tively, it can be considered as a scalar field defined 
on bonds). The function 

1 1 

F ^ X > ~ 2x\n(J (x)/J 1 (x)) 2x\n{A/xY X < ' 

is introduced to "correct" the Villain transforma- 
tion for small Ej (see e.g. Ref. Q). As we see, 
its entire effect is to renormalize (increase) the 
Josephson coupling Ej — > EjF(eEj), but it does 
not affect the physics. In the following we will 
use only the renormalized constant. 

The rhs. of Eq. (Al) can be rewritten as 



— iJj 



rVfar} (A2) 



Now the Gaussian integration over the phases can 
be easily performed, yielding 

z = exp { ~ 2e2e ]l i^v 1 ^ 

i,r 

and the summation is constrained by the conti- 
nuity equation, 

VJ lT - q iT = . 
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The time derivative stands for a discrete deriva- 
tive /(r) = e- x [/(r + e^) - /(r)]. The constraint 
is satisfied by the parameterization [Q 

Here the operator (nV)" 1 is the line integral 
on the lattice (in Fourier space it has the form 
i(k x + fey) -1 ), e^"' is the antisymmetric ten- 
sor, while Ai T is an unconstrained integer-valued 
scalar field. 

With the use of the Poisson resummation 
(which requires introducing a new integer scalar 
field Vi T ) the partition function can be rewritten 
as 

Z = ^2 exp-S{q,v} . 

qiT, ViT 

The effective action for the integer charges qi and 
vorticities Vi is 

S{q,v}=\2e 2 eT l i j rQirOr j 1 q j r 



2eEj 

eE 3 ^ r _ 



— J2\r 



4tt 



ijT 



2nvi 



eE 3 



^(nV) -1 ©. 



Gi 



eEj 



(A4) 

The kernel Gy is the lattice Green's function, i.e. 
the Fourier transform of k~ 2 . Finally, after some 
algebra |2?J we arrive to the effective action (p"2|), 
rewritten in the continuous notations. 
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